Fermion Cooper Pairing with Unequal Masses: Standard Field Theory Approach 
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The fermion Cooper pairing witli unequal masses is investigated in a standard field theory ap- 
proach. We derived the superfluid density and Meissner mass squared of the U{1) gauge field in 
a general two species model and found that the often used proportional relation between the two 
quantities is broken down when the fermion masses are unequal. In weak coupling region, the su- 
perfluid density is always negative but the Meissner mass squared becomes mostly positive when 
the mass ratio between the pairing fermions is large enough. We established a proper momentum 
configuration of the LOFF pairing with unequal masses and showed that the LOFF state is energet- 
ically favored due to the negative superfluid density. The single plane wave LOFF state is physically 
equivalent to an anisotropic state with a spontaneously generated superflow. The extension to finite 
range interaction is briefly discussed. 
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I. INTRODUCTION 

The asymmetric cooper pairing between different 
species of fermions with mismatched Fermi surfaces, 
which was discussed many years ago, promoted new inter- 
est in both theoretic and experimental studies in recent 
years. The mismatched Fermi surfaces can be reahzed, 
for instance, in a superconductor with Zeeman splitting 
induced by an external field[l|,Qiyil3i an atomic fermion 
gas composed of two species of atoms with different den- 
sities and/or masses [algi, an isospin asymmetric nuclear 
matter with proton-neutron pairing 0, and color super- 
conducting quark matter with charge neutrality[a,|9|,|l0|. 
Among the mechanisms which can produce asymmetry 
between the pairing fermions, the mass difference is a 
very robust one. The cooper pairing between fermions 
with unequal masses was firstly investigated by V.Liu and 
F.Wilczek|l2|. They considered a fermion gas composed 
of light and heavy fermions with attractive interaction. 
A homogeneous and isotropic pairing state which is sim- 
ilar to the Sarma state [j] was proposed to be the ground 
state of such systems. Such an exotic pairing state is now 
called breached pairing (BP) state or interior gap state. 
In the BP state there exists gapless fermion excitations, 
and the superfluid Fermi gas and the normal Fermi gas 
coexist in the momentum space. 

It was found many years ago that the Sarma state 
suffers a thermodynamic instability 1]. It is now gener- 
ally accepted that this Sarma instability can be cured 
in some physical conditions, such as a long-range in- 
teraction when charge neutrality is required [3, Ii3j a 
proper finite range interaction between the two species 
of fermions with large mass difference|l2j|. and a super- 
fluid Fermi gas with density imbalance in strong cou- 
pling region jlSl Il4l llSl ll q | . While the Sarma instability 
can be cured, it was soon found that the superfluid den- 
sity of the BP state is negative 17] and the free energy 
of the mixed phase is also lower than that of the BP 
state 18, 19]. Meanwhile, in the study of color supercon- 
ductivity, it was found that the gapless color supercon- 
ductors possess paramagnetic response to external color 



magnetic fields, i.e.. the Meissner masses of some gluons 
are negative|2fl El EHH Hi]. AU these phenomena 
indicate that the homogeneous and isotropic BP state is 
unstable and some spatially inhomogeneous and isotropic 
states are energetically favoredj2ai2al27Ll2ll2ll30ll3ll]. 
The superfluid density is a fundamental quantity in su- 
perconductivity. It is well known that the superfluid den- 
sity is proportional to the Meissner mass squared and 
often measured via the London penetration depth in ex- 
periments. Due to this relation, one often regards the 
negative superfluid density observed in a BP state and 
the negative Meissner mass squared observed in a gapless 
color superconductor as the same instability. 

However, almost all these studies focus on systems 
where the masses of the pairing fermions are equal. In 
the study of superfluid stability of interior gap states, 
S.T.Wu and S.Yip derived a formula of the superfluid 
density for non-relativistic asymmetric fermion superflu- 
ids with the concept of a uasiparticles \ll\ . In the equal 
mass case, their formula is consistent with the result cal- 
culated from the linear response theory [32], the current- 
current correlation functionl33j, and the Meissner mass 
squared 34]. However, in unequal mass systems, the for- 
mula of superfluid density seems quite different from the 
Meissner mass squared JS^. Does the proportional re- 
lation between the superfluid density and the Meissner 
mass squared still hold in unequal mass systems? In 
this paper, we will derive the superfluid density and the 
Meissner mass squared in unequal mass systems in a stan- 
dard held theory approach where the superfluid density 
and the Meissner mass squared are treated in the same 
way. 

The paper is organized as follows. In Section ^ we 
briefly review the formalism of the two species model. In 
Section IIIII we derive the formula of the superfluid den- 
sity and compare it with S.T.Wu and S. Yip's phenomeno- 
logical method. In Section IIVI we derive the Meissner 
mass squared and show that it is not proportional to the 
superfluid density in unequal mass systems. The super- 
fluid density and Meissner mass squared in the breached 
pairing states are calculated in Section The LOFF 



pairing in unequal mass systems is discussed in Section 
IVil The extension to finite range pairing interaction is 
briefly discussed in Section IVIII We summarize in Sec- 
tion IVIIII We use the natural unit of c = ft = fc^ = 1 
through the paper. 



II. 



TWO SPECIES MODEL 



The physical system we are interested in in this paper 
is an idea system composed of two species of fermions 
with attractive interaction. The system is described by 
the Lagrangian density with imaginary time r = it, 

^=Y1 ^'' ( "^^ + ^~ + M ^* + 9^*a^t^bi^a, (1) 



i^a^b 



where ipi = ipi (x) with x — (r, x) are fermion fields for the 
two species a and 6, the coupling constant g is positive 
to keep the interaction attractive, ma and m?, are the 
masses for the two species, and fj,a and /if, the chemical 
potentials. 

The key quantity to describe a thermodynamic system 
is the partition function Z. It can be expressed as 



[di;,][dr^Wod-Id'^^li'^,4':] 



(2) 



in the imaginary time formalism of finite temperature 
field theory, where P is the inverse of the temperature T, 
(3 = 1/T. For attractive interaction g, we can perform an 
exact Hubbard-Stratonovich transformation to introduce 
the auxiliary boson field (t>{x) and its complex conjugate 
(j)*{x). With the Nambu-Gorkov fields ^, 'I' defined as 



^(a;) 



tpa 



^{X) ^{ra i'b) , 



(3) 



we can express the partition function as 

[d^] [d^ [d^] [d^*]e^o dr I d^x(*K;*-|0| Vs) ^ (4) 



where the kernel /C[(/), (f>*] is defined as 



-9r + ^+^^a 



-9r-^~ Mb 



(5) 



In mean field approximation, we replace and (/>* by 
their ensemble averages A and A*. In a homogenous 
and isotropic state, they are independent of coordinates. 
Then we can directly evaluate the Gaussian path integral 
and obtain the thermodynamic potential 



n 



\A\- 
9 



^E/(^Trlng-i(zc.„,p) (6) 



in terms of the inverse fermion propagator 

V A' lUJn + f^ 



(J) 



with the free fermion dispersions cL = p^/(2mi) — /i^. 
The explicit form of the fermion propagator which we 
need in the following sections can be explicitly expressed 
as 

^(^c.„,p)=f^";^""'Pj ^i^S^'^-'P)) (8) 



iLU„ - eA + fs 



with the matrix elements 

5ii(icj„,p) = 

^22(JWn,p) = 



52l(JWn,P) = J-. To T' (9) 

[UVn - e-A) - ^A 

where the quantities eg, e^ and ea are defined as 



{iujn ' 

iuJn 


-^aY 


-4' 


{iuJn - 


-^aY 

-A 


-4' 


{iLU„ - 


-A* 


-4' 



es,A^{el±el)/2, e^ = Jel + \A\ 



(10) 



Since all quantities depend only on |A|, we can set A to 
be real from now on. From the pole of the fermion propa- 
gator we can read the dispersions Cp and e^ of fermionic 
quasiparticles: 



.A 



"P 



(11) 



For eA = 0, we recover the well know BCS type exci- 
tation. The asymmetric part sa is the key quantity to 
produce exotic pairing states. 

The occupation numbers of the two species of fermions 
can be calculated via the diagonal elements of the fermion 
propagator. 






Uaip) = Tlini > Cii(zw„,p)e 



Hb 



n 

(p) = -riim V'c;22(icj„,p)e" 



(12) 



Completing the Matsubara frequency summation, we ob- 
tain 



na(p) = ulf{e^)+vlf{~el), 



(13) 



with the coherent coefficients vl = (l + es/eA)/2 and 
Wp = (1 — es/eA)/2. The particle number densities n^ 
and nf) for the species a and b are obtained by integrating 
na(p) and ^^(p) over momentum. 



III. SUPERFLUID DENSITY 

In this section we try to derive the superfluid den- 
sity in a standard field theory approach. When the su- 
perfluid moves with a uniform but small velocity Vs, 



the condensates transform as A ^ Ae^'^'^'^ and A* -^ 
^*g-2iqx ^[f^iy ^j^g total momentum of the cooper pah- 
2q = (jria + TOfc)vs, and the fermion fields transform as 
■0a -^ ipae^"^"^'^ and i/'b -^ tpbC^'^'''^ with the momenta of 
the two species Qa — rUa'Vs and qb = mb'Vs which satisfy 
Qa + qh = 2q. The superfluid density tensor pij is defined 
as 113 



n{vs) ^ n{o) + j. • V, + -Ay (v,),(v,)j- + • • • . (14) 



For a homogeneous and isotropic superfiuid, we have 
Pij = (5ijPs/3, and the above formula can be reduced 
to 



n(v,)=f)(0)+j, -V, + -p,v2 



(15) 



where p^ is the superfluid density. When ps is nega- 
tive, the homogeneous and isotropic state is unstable and 
a state with spontaneously generated superflow which 
breaks the rotational symmetry is energetically favored. 
After the transformation of the condensates and the 
fermion fields, the thermodynamic potential is changed 
as 



f](v,) = — -T^ /'^Trlng7i(*c.„,p 

in terms of the Vj-dependent inverse propagator 

A 



) (16) 



VjJn - e 



Qs («^n,P) - I ^ 

Using the relation 



iujn + e' 



p-qt 



(17) 



Q's = ^^ - P • V, - -i]„ 



(18) 



with the matrix E.^ = diag(ma,~rni,), we can do the 
derivative expansion 



Tr In g-^ - Tr In g-i = p • v,Tr {G) - ^Tr (gS„0 



-(p-v,)2Tr(5g) + ....(19) 



With this expansion, we can expand r2(vs) in powers 
of Vs. The superfluid density can be read from the 
quadratic term in v^. After some direct algebras, we 
obtain 



Ps = maUa + mbUjj 



(PTP P^ 



(rn+r22 + 2Ti2)(20) 



(27r)3 3 
where 7ii,722,7i2 are the fermion Matsubara frequency 



summations deflned as 



9 2 
UpVp 



f{4) + f{e^)-l 



EA 



(21) 



+Kr{4)+vin4), 



%2 = Ty ^^22^2 



2 2 

UpVp 



.fi4) + f(4)-^ 



4 j//,B\ 



EA 



4f'i4)+4n4 



'?12 = r / ^^12^21 



2 2 



EA 



/'«) + /'«) 



with f{x) being the Fermi distribution function and 
f'{x) = df{x)/dx. Using these results we get 



maTla + TObTlfc + 



d3 



p p 



(27r)3 3 



[f\e^) + f'{4)].{22) 



One can easily check that this formula is invariant under 
the exchange a ^^ b. When A = 0, i.e., in the normal 
state, we have 



From the identity 

/■oo „4 i-oo 2 



(23) 



(24) 



Ps vanishes automatically in the normal state. 

The result we obtained here is in agreement with the 
formula derived by S.T.Wu and S.Yip with a phenomeno- 
logical method[l7|. With their method, in presence of a 
small superfluid velocity v^, the quasiparticle energies 
are shifted by p • Vg and the occupation numbers become 

fiaiv) = Mp/(ep+p-v,)+t;2/(-e^ + p.v,), 
n6(p) - ulf(e^+p-ws)+vlf(-e^ + p-ws).{2?,) 

The number current can be decomposed into a diamag- 
netic and a paramagnetic parts: 

d^V ~ , . d 

rn,(p)v,. = Pi^s, 



(27r)3 



m,. 



J2n)3 



pn,(p) = pfvs. 



(26) 



To leading order in Vs, we have pf = Ui. Using the fact 
J? = for Vs = 0, i.e., 






d^p 

(27r)3 



p(ni(p)-n,(p)). 



(27) 



we obtain 



Pb = — 
nrib 






(27r)3 3 



bin4)+4n4)\- (28) 



The total superfluid density is defined as ps — niaPa + 
rUhPh with pi = pf + p^^. Using the fact u^ + Vp = 1, 
it is exactly the formula we obtained above. From our 
derivation, we can also decompose the superfluid density 
into two parts 



p^ = maUa + 



Pg = nibfib + 



l2^"3 



d3 



P P 



(27r)3 3 



(rii+Ti2), 



(722 + 7i2) 



(29) 



where p° = rUaPa and p\ = rribpi, can be defined as the 
superfluid densities for the two species of fermions. 

One may ask why these two derivations give the same 
result. In presence of a small superflow Vs, the quasipar- 
ticle dispersions can be read from dett/~^ — 0. After a 
simple algebra, we get the modified dispersions for the 
quasiparticles 



4 



A2 



?2 



CA+P-^S 



CA - P ■ Vs 



(30) 



with 



J- , -2 

es = es + jKi^a + 'mb)v^, 
eA^ eA + -7 {ma - rnb)wl. 



(31) 



To leading order in v^, the quasiparticle energies are re- 
ally shifted by p • v^. However, one should note that the 
derivation with the concept of quasiparticle is inconsis- 
tent. The particle occupation numbers in presence of a 
superflow should be 

na{p) = u^/(e^ + p-v,) + i;2/(-e^ + p.v,), 
nb{p) = ulf{e^-p-v,) + vlf{-e^-p-v,) (32) 

to leading order in Vg. Using this correct occupation 
numbers, one can not obtain the correct result since p^ 
will change a sign. We guess that for such an asymmetric 
system, one can not self-consistently derive the superfluid 
density with the concept of quasiparticles. Only for sym- 
metric systems like the standard BCS, the method works 
as discussed in text books. 

The formula of the superfluid density we derived here 
is in principle only suitable for grand canonical ensembles 
where the chemical potentials p,a and pb are fixed|38|. For 
systems where the particle numbers n^ and Ub axe fixed 
or the total number n ~ Ua+nb is fixed, we should use the 
free energy T = il + p,ana + Pbnb or JF = fl + (pa+Pb)'n/2 



instead of the thermodynamic potential fi to calculate 
the superfluid density. However, we can show that such 
a correction is beyond order 0{v'^) for an homogeneous 
and isotropic state[39j, and hence we can safely apply the 
above formula to the systems with flxed particle numbers. 
For instance, if n^ and Ub is fixed, we have 



52^ 



9v2 






dn.i dp,i 



(33) 



v.,=0 



For a homogeneous and isotropic state, the term 
i9ni/i9vs|vj=o vanishes automatically. 

In our derivation, we did not use the assumption of 
weak coupling, and the formula can be applied to study 
the superfluid stability in the BCS-BEC crossover in a 
light-heavy fermion gas such as a mixture of ^Li and ^"K. 
In recent studies on BCS-BEC crossover in equal mass 
systems, it was found that the BP state is stable in the 
BEC region, i.e., it is free from the Sarma instability 
and negative superfluid density[l3. We expect that such 
a stable BP state can also be realized in a light-heavy 
fermion gas in strong coupling. 



IV. MEISSNER MASS 

The two species model is invariant under the following 
phase transformation 

tpi(x) -^ e'^^^^j^ix), (j){x) -^ e*(^"+'^''V(a;) (34) 

with arbitrary and constant phases (pa and (pb, which 
means that the symmetry group of the model is C/(l)<p„ <Xi 
U{1)^^. The order parameter is invariant only for ipa — 
—ipb- In presence of a nonzero expectation value of 0, 
the symmetry group is spontaneously broken down to a 
U{1) subgroup 



U{l)^^^U{\)^,-^U{\) 



Va-Vb- 



(35) 



The unbroken f/(l) subgroup corresponds to the phase 
difference A(/5 = i^a — 'Pb, and a Goldstone mode corre- 
sponding to the total phase ip = ipa + ^b will appear. 
Let's add a U{1) gauge field A^ in the Lagrangian, 



i^a.b 



E^U-Dr. + ^rr 



2mi 



with Dfj,i = 9^ - leQiA^ 



Pi \ -ipi + gt/j^ipl^Jbi^a + Ca 

(36) 
where La is the gauge field 
sector, and eQa^ &Qb are the gauge couplings for the two 
species of fermions. In presence of a gauge field, the Gold- 
stone mode disappears and the gauge field will obtain a 
mass TTiA via Higgs mechanism. This is called Meissner 
effect in superconductivity, and the mass ttia the gauge 
field obtains is called Meissner mass. 

The standard way to calculate the Meissner mass is 
to evaluate the polarization tensor H^^ iK) of the gauge 
field. For the interaction (|36|l . the spatial components of 
the polarization tensor read 



n,,(fc) = Hf,(fc) + Hf,(fc) 



(37) 



with the diamagnetic part 



,T 



nf,(fc) = <J.,e2-^Tr[g(p)Erf] 



and the paramagnetic part 



n?,(fc) 



V 



(38) 



^ PiPjTr [t; (p+) Sp^ (p_ ) Ep] (39) 



The coefBcients Mf, are just the components of the 
Meissner mass squared tensor. In a homogenous and 
isotropic superconductor, we have Mf^ = for z 7^ j 
and Mil = -^22 = -^33 j ^^^ the Meissner mass squared 



m'ji is defined as 



^Ml 



(44) 



with p± — p ± fc/2, where the matrices T,d and Sp are 
defined as 



The thermodynamic potential in presence of the static 
and long wave vector potential A can be expressed as 



91 

rria 

-S 




(40) 



Note that the paramagnetic part is just the current- 
current correlation function which gives both dia- 
magneti c p art and paramagnetic part in relativistic 
systems I20I l2l| but only paramagnetic part in non- 
relativistic systems. The Meissner mass mA can be eval- 
uated by 



m\ = - lim((5y - fcifcj)!!^ (0,k). 



(41) 



If m^ is negative, the homogeneous and isotropic state 
suffers the magnetic instability |23, UM ^^^ ^ state with 
gauge field condensation (A) ^ which breaks the rota- 
tional symmetry is energetically favored. 

For a clear comparison of the superfluid density and the 
Meissner mass squared in unequal mass systems, we em- 
ploy another approach 34j. The Meissner mass squared 
can be calculated via the response of the effective poten- 
tial to an external transverse vector potential. In pres- 
ence of a small external vector potential A(0, q ^ 0) in 
the static and long wave limit, the effective potential of 
the system can be expanded in powers of A, 



n{A) = n(o) + 3a-a + -m^a ■ a 



with the coefficients 



Mf. 



d^n{A) 



dASA, 



A=0 



(42) 



(43) 



Q.{A) 



A2 



^E 



dh 



9 —J (27r)3 

with the A-dependent inverse propagator 



Trln^^i(iw„,p) (45) 



''p-eQbA 



(46) 



Using the same trick in Section ITTll we have the relation 

(47) 
and the derivative expansion 



^^i = g-i-eSpp-A--ErfA2 



Tr In g^i - Tr In ^^^ == ep • ATr (^Ep) 
-yA^Tr (gSrf) - y(p • A)2Tr (^Ep^Sp) 



(48) 



The Meissner mass squared can be read from the 
quadratic terms in A. After some algebras we obtain 



m\ 



m„, " nib 
2 / d'p p^ (Ql, 

^2 



^'^—Ql + —Ql 



e- I -^-^-r I ^Tii + ^T22 



(49) 



(27r) 



m" 



2QaQt 

marub 



Ti: 



The second term is just the long-wave and static limit of 
the current-current correlation function. With the result 
of frequency summations in Section IIIII we obtain an 
explicit expression, 



m\ 



m„ 



mb 



-Ql 



Wb 
mb 



UpVp 



d^p p^ 

(27r)3 7 



2..2/(^p)+/(^g)-l 



fA 



Qa 



Qb 



mb P 



/'«) + 



The formula is invariant under the exchange a ^^ h. For A = 0, it is reduced to 



Qa 



Qb 2 

— S 



no 



(50) 



dp 



27r2 



rua ° mb 



dp 



67r2 



f'i4). 



+ 



i_l_Pir)2 

~~~2~'^b 



(51) 



and vanishes in the normal state, as we can expect. 

The formulae of the superfluid density and the Meiss- 
ner mass squared we derived seem quite different. In the 
symmetric case with rua — rni, = m and Qa = Qb = 1, 
we recover the well know result 133, ^3 



Ps 



mn - 



d^ 



(27r)3 3 



^ m m^ J (27r)3 3 

and the proportional relation 

Ps m"^ 



[/'«) + /'(e^)],(52) 
(53) 



at any temperature T < Tc- In fact, one can easily check 
that for systems with Qajvna = Qh/wh, the proportional 
relation still holds. However, for general asymmetric sys- 
tems, this relation is broken down. 



V. THE BREACHED PAIRING STATE 

We calculate the superfluid density and the Meissner 
mass squared for the breached pairing state at zero tem- 
perature in this section. Explicitly, the dispersions of the 
fermionic quasiparticles can be expressed as 




A^±(£^ + ^^ 



(54) 



with the reduced masses m = 1mam\,j(rno, + mb) and 
m! = 2mami, / {rrifj — ma) and chemical potentials ^ = 
{Pb 



{Pa + Pb)/'2 and Sfi 
that for 



/ia)/2. One can easily check 



A< Ac 



\pI 



pI\ 



iy/niamb 



(55) 



with Pi = ^y2mifii, one branch of the fermionic quasipar- 
ticles can cross the momentum axis and hence becomes 
gapless. This is the so called breached pairing state or 
interior gap state |ll|. The gapless nodes determined by 
Ep = or Ep = are located at p = pi_2 with 

1 



Pi, 2 — 



pI 



(pI ~ Pb)'^ - 16maTOhA2. (56) 



The gap equation which determines A at zero tempera- 
ture reads 



Anas 



(27r)3 



6(6^) - 9(-6^ 

2eA 



1 



(57) 

Here Q{x) is the step function, and the s-wave scattering 
length Us is related to the bare coupling g via 



1 



47ra, 



+ 



d^p 1 



(58) 



with ^p = p^ /2m. One can employ another regulariza- 
tion scheme where the pairing interaction is restricted 
in a narrow region around the common Fermi surface 
Pf — ^/2mjl 11] which is suitable only for weak coupling 
case. We have checked that the results from these two 
regularization schemes are the same in weak coupling. 

In the BCS phase, all fermionic excitations are gapped, 
the gap equation can be reduced to 



ATTU:, 



dp 



1 



^(2 

27r2 l^2eA 2^ 



and the solution in weak coupling is 



Ao 



Be ^/ie 2pFio 



(59) 



(60) 



Since all fermionic quasiparticles are gapped, we have 
na — TT-f, — n/2, and the superfluid density reads 

Ps ^ maHa + mbUb = -{ma + mb)n, (61) 

which means that all fermions participate in the super- 
fluid. In weak coupling, we have approximately n ~ ^^ 
and 



d^p p2 ^2 mp^ 



(27r)3 12 e| 127r2 



which leads to 



2 ne^iQa + QbY 
m^ - 



2{ma + mb) 



(62) 



(63) 



Note that m^ depends only on Qa + Qb, and this con- 
clusion is valid only in weak coupling we considered by 
hand. 

In the breached pairing state with gapless fermionic 
excitations, the gap equation can be expressed as 

m_ _ r , p^ (J- \_\ _ r^ , p' 1 

~47Ta,^ Jo ^ 2^2 (^ 2e A " 2ep j " ip, * 27r2 2e a ' 

(64) 
Using the equation for the BCS gap Aq, the solution of 
the gap equation can be well described bv,19.| 



A = VAo(2A, -Ao). 



(65) 



The gap varies in the region < A < Aq, and the differ- 
ence between the Fermi momenta satisfies 



2y/m^mbAo < \Pa ~ Pb\ < 4:^/rrwrnbAo. (66) 

If we define the density asymmetry a = \na — nb\/{na + 
Ub), then A = Aq corresponds to a = and A = cor- 
responds to the maximal asymmetry ac- When a varies 
from to ttc, the gap A decreases from Aq to 0. 

We firstly discuss the case rua — m,b where the su- 
perfluid density is proportional to the Meissner mass 
squared. In this case we have Ac = Sp, and p'^ /2m' = 0. 



Let's set Sfj, > without loss of generality. At zero tem- 
perature the superfluid density reads 



Ps = rrm 



1 

6^ 



dpp'^S{eA - Sfi). 



(67) 



In the weak coupling region, eA — Sn — has two possible 
roots pi , p2 and ps can be evaluated as 



Ps = mn 



1-77 



5iiQ{5p - A) 



(68) 



with the coefficient 77 = (pf + P2) / (67r^?T,) . Since the 
coefficient 77 is approximately equal to 1, ps can be well 
approximated by 



5^iQ{5p - A) 



A2 



(69) 



It is now clear that in the BP state with A < 5p, ps be- 
comes negative and the BP state is unstable. We should 
emphasis that the function in the bracket is universal for 
gapless fermion superfluids in equal mass systems. This 
function appears in the Meissner mass squared for the 8th 
gluon in two flavor gapless color superconductor ^20. .211] . 
In some anisotropic states in equal mass systems such 
as the LOFF state[2a, |33 and the BP state via p-wave 
pairing|43, a similar function appears. In the LOFF 
state ^u is replaced by an angle dependent mismatch 
b^ [2^ I33 , and in the p-wave BP state the gap A is re- 
placed by an anisotropic gap function An 40] . 

In general case with ma 7^ f^Zb, we define a mass ratio 
A — mb/rria and set A > 1 without loss of generality. For 
X ^ 1, the results for pa < Pb and Pa > Pbior ria < Ub 
and Ua > Ub) are not symmetric. We will discuss these 
two cases separately at zero temperature. For the sake 
of simplicity, we set Qa = Qb = 1 in the calculations. 



A. Pa < Pb{na < rii) 

In this case, the branch e^ becomes gapless and we 
have Uaip) — 0, ^^(p) = 1 in the region pi < p < p2- At 
zero temperature, the superfluid density in the BP state 
can be evaluated as 



Ps = rua 



67r2 



(70) 



with the coefficients Us and [3s defined as 
A 



ii-(A+i).?r'^'+^)/./^?f^- 



Ps = \- 



+ 3(A + 1) / dp^vlijl) 

JR-2 P2 



\1-{X + I)vl\ ' "'- ' ^M„ ">^.3«P 



where Vi and i"! are the values of vi sip = px and p — P2, 
and the integral regions i?i and R2 are < p < pi and 



P2 < p < 00 respectively. The Meissner mass squared in 
the BP state can be evaluated as 



mb OTT'^ 

where the coefficients a^m Pm and 7„i are defined as 



(72) 



[l + {\-l)viY f ,P' 2 

"^ = ^- \l-iX + l)vf\ -3(A + i)y^^dp^v 



Pm = 1 



[1 + {X-I)vj 



|1-(A + 



>|.3(A.l)/^^./-!.^, 



(A -1)2 



dp 



A% 



A + 1 Jr,+r, V^[(ep-M)^+A2]3/2 



,(73) 



From Pa < Pb, we have Xpb > Pa, and the chemical 
potentials in the BP state satisfy 



Ao Xpb- Pa ^ . 



(74) 



Without loss of generality, we can keep pt fixed. After a 
simple algebra we find that 

A^b - 2\/AAo < Pa < Xpb - VAAq. (75) 

The lower bound corresponds to A = Aq where a — 0, 
and the upper bound corresponds to A = where a = 
ac- Then we can calculate the superfluid density and the 
Meissner mass squared as functions of A/Aq in the BP 
range < A/Aq < 1. In FigQ] we show the superfluid 
density and Meissner mass squared for different values 
of mass ratio A. We found that the superfluid density 
is always negative at any mass ratio, but the Meissner 
mass squared is positive in the region < A/Aq < v 
with V < 1. When the mass ratio becomes very large, 
such as A = 100, v is close to 1. Even though there 
exists a big room where the Meissner mass squared is 
positive, both the superfluid density and the Meissner 
mass squared tend to negative infinity near A/Aq = 1. 
Such a divergence at the BP-BCS transition point, which 
comes from the divergent density of state of the gapless 
excitations, can not be avoided [l7ll34l |. 



B. Pa > Pb{na > nt) 

In this case, the branch e^ becomes gapless and we 
have ria(p) — 1, nb{p) = in the region pi < p < p2. At 
zero temperature, the superfluid density in the BP state 
takes the same form (|70|l but with different coefficients 



as ^ X- 

Ps = l- 



7)2 
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FIG. 1: The scaled Meissner mass squared Fa ~ 
m,^ /mo (solid lines) and superfluid density Fs = ps/po (dashed 
lines) as functions of A/Ao for different values of mass ratio 
A in the case Pa < pt- The normalization constants mo and 
po are chosen to be ttiq = e m Pf,/{mamb) and pa = rribPt- 
The BCS gap Ao is chosen to be Ao = O.Olpb- 



The Meissner mass squared takes also the same form (|72|l , 
but the coefficients «,„, Pm are modified to 



= 1 



[X-{X-l)vl 



|A-(A + 



^^-3(A+i)y^^dp-3.,, 



}vi] 



|A-(A + l)i'^| Jr^ p^ p 



P\.2 



(77) 



and 7m remains unchanged. 



FIG. 2: The scaled Meissner mass squared i^A (solid lines) and 
superfluid density Fs (dashed lines) as functions of A/Ao for 
different values of mass ratio A in the case pa > Pb- The BCS 
gap Ao is chosen to be Ao = O.Ol/ib. 



Similarly, for fixed /ib, we have 

X^ib + VAAo < Ha< Xnb + 2\/AAo (78) 

for the BP state. The superfluid density and the Meissner 
mass squared are calculated in Fig|21as functions of A/Ao 
in the range < A/Aq < 1. The qualitative behavior is 
almost the same as in the case pa < Pb- 

In summary, we have shown in weak coupling that, 
the superfluid density of the BP state is always negative, 
while the Meissner mass squared can be positive in a wide 
region. The conclusion here is valid for stronger coupling, 
if there exist two gapless nodesjla|. 



VI. LOFF STATE 

In weak coupling the superfluid density of BP state 
is always negative, which indicates that the BP state is 
unstable and some inhoniogeneous and anisotropic state 
is energetically favored. In this section we show that the 
LOFF state is energetically favored due to the negative 
superfluid density. 

For the sake of simplicity, we consider the simplest 
pattern of LOFF state, namely the single plane wave 
ansatz 



i^x)) ^ Ae'^^-- , (r(x))=Ae- 



2«q-x 



(79) 



where A is a real quantity, and 2q is the so called 
LOFF momentum which is the total momentum of a 
Cooper pair. To evaluate the thermodynamic potential 
of the LOFF state, we often define new fermion fields 
Xaix) = e'^'^'^tpaix) and Xbix) — e^'^'^^pbix), and then 
can directly evaluate the Gaussian path integral in the 
new basis XayXh- Following this way, the thermodynamic 
potential reads 



A2 

9 V 






TTlng-\iLUn,p) 



(80) 



in terms of the new inverse propagator 



g \iujn,p) 



A 



p+q 



A 



This is just the (q + p, q — p) picture of the LOFF pair- 
ing, which means the fermions in the cooper pair move 
together with a total momentum 2q. To see whether the 
LOFF state is energetically favored, we take the small q 
expansion 



1 (9^0 



q' + Oiq% 



(82) 



q=0 



where we have chosen a suitable z-direction such that 
q = (0, 0, q). Notice that a linear term in q vanishes au- 
tomatically. One can easily observe the following relation 
between the momentum susceptibility and the Meissncr 
mass squared: 



,d^n 



q=0 



(83) 



For large mass difference, since the Meissner mass 
squared is positive in the small A region which is just 
the window for LOFF state, as we have shown in the 
last section, we may conclude that the LOFF state is 
not energetically favored. However, we shall argue in the 
following that this is not the truth. 

To give a correct argument we focus on the fact that 
the superfluid density is always negative for any mass 
ratio. Notice that we can do any transformation like 

Xaix) = e'^^-^M^), Xbix) ^ e'^'^-^Mx) (84) 



to evaluate the effective potential, since the phase factor 
in the condensate can be eliminated by any q^ and qf, 
satisfying qa + qt = 2q. For such a general transforma- 
tion, the thermodynamic potential reads 



A2 

n^ — 

9 



T 






dh 



(27r 



-Trlng^_^^^^(iw„,p) 



with 



■^p+qa 



p-qb 



(85) 



(86) 



This arbitrariness of phase transformation is directly 
linked to the fact that the symmetry group of the model 
Lagrangian is U{1)^^ ® U{1)^^. To link the LOFF state 
and the superfluid density, we introduce a LOFF velocity 
w such that 



qa = m^w, q;, = mfcW. 



(87) 



With a suitable choice of coordinates such that w = 
(0, 0, w), we can do the similar small w expansion 



n{w) - n{Q) ^ - 



id'^n 



2dw^ 



+ 0{w^). (88) 



Also, one can easily observe the following relation be- 
tween the velocity susceptibility and the superfluid den- 
sity: 



d^n 



Ps 



dw"^ 



(89) 



This intuitive argument indicates that the energetically 
favored momentum configuration of LOFF state is (\a = 
ma'w, qf, = nib'w, which means the single plane wave 
LOFF ansatz is nothing but the anisotropic state with 
spontaneously generated superflow v^ if we identify w = 
Vg. In fact, it is quite easy for us to understand this 
fact. The physical picture of the LOFF state is that the 
fermions in a cooper pair move together with a nonzero 
momentum, and hence they should possess a same ve- 
locity, not momentum. In fact, we have checked nu- 
merically that for general choice of q^ and q;, that the 
quadratic term in the expansion is always negative only 
when qo = moW, qf, = mbVf. 

With the proper configuration qa = rna'w, qt, = nib'w, 
we can evaluate the effective potential as 



n{A,w) 



A2 



/(^^^^-^^) 



/(^[^(^p) + ^(^p)]- (90) 



Here H{x) = Tln{l + e"^/^), Es = ^p - fJ-n,, Ea 



y^E'g + A^ and E^, E^ are the energies of the quasipar- 



ticles 



P 



e-^- = Ea±[^ + s^ + p-^ 



2m' 



(91) 
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with fij^ = fi — {ma + mb)w'^ /A and Sw = Sfj, + {rUa — 
mt)w^/4:. Notice that a new kinetic energy term {nia — 
mb)w'^/4: which vanishes in the equal mass case appears 
in the asymmetric part of the quasiparticle dispersions. 
In fact, this momentum configuration is the most con- 
venient one to calculate the LOFF solution, since the 
anisotropic term p • w appears only in the asymmetric 
part. 

Now we give a preliminary discussion on the LOFF 
state with mass difference. For convenience, we assume 
the pairing interaction is restricted in the region pp — 
A < \p\ < pp + A with A ^ pp, here A serves as a 
natural ultraviolet cutoff in the theory. In weak coupling 
we can safely neglect the terms of order 0{w^) and do 
the following replacement 



p • w ^ ppw cos 9, 



P 
2m' 



2m' 



(92) 



where 9 is angle between p and w. Up to now, all things 
become the same as those in the equal mass systems j^ 
l2q| . and the conclusions there can be directly applied. If 
the chemical potentials for the two species are fixed, the 
corresponding LOFF window is 



0.707Ao < 



2m' 



< 0.754Ao, 



(93) 



where Aq is the BCS gap , an d the LOFF velocity w is 
approximately given by[j,|26j 



ppw ~ 1.2 



X ^ Pf 



(94) 



Defining the mass asymmetry e = {mi, — ma)/ {mb + ma) 
we have 



5/i - 



2m' 



1 



[(l + e)Ai6-(l-e)/ia] = 5(e), (95) 



and then can reexpress the LOFF window as the conven- 
tional form in equal mass case 



0.707Ao < |(5(e)| < 0.754Ao 
and ppw ~ 1.2(5(e). Due to the relation 



(96) 



5{e) 



9 ? 

Pb-Pa 



2{ma + mb)' 
the size of the LOFF momentum is 

|2q| = {ma + mb)w ~ l.2{pb ~ Pa), 



which is just we expect. A LOFF state induced by a pure 
mass difference is of great interest since in some physical 
systems the chemical potentials are always equal due to 
chemical equilibrium. Setting ^^ = /Xf, = /i we obtain the 
mass difference window 



0.707— < lei < 0.754—. 



(99) 



In weak coupling, Aq <C /x, the LOFF state can exist 
only when the mass asymmetry is very small, which may 
be realized in electronic systems. If the particle number 
densities tIq and rif, are fixed, we should compare the free 
energy J^ = Q + ^a^a + fJ-b^b- In this case we have the 
similar expansion !F{w) — T{0) + psW^ /2 + 0{w'^) which 
means the LOFF state is more stable than BP state. The 
LOFF window will be larger, which is similar to the equal 
mass system[3^|33- Such a situation may be realized in 
cold atomic Fermi gas, such as a mixture of ^Li and '^"K 
atoms. 

Finally, we calculate the superfluid density tensor and 
the Meissner mass squared tensor in the LOFF state. 
Since the rotational symmetry 0(3) is broken down to 
0(2), the superfluid density and Meissner mass squared 
become tensors pij and {m\)ij. We can decompose them 
into a transverse part and a longitudinal part 

Pi] = Pt{Sij - WiWj) + PLWiWj, 

{m\)ij = m^{6ij - WiWj) -f m\wiW.j (100) 

with w = w/|w|. The transverse and longitudinal super- 
fluid density read 



PT = maUa + mbUb + 



PL = mana + mbUb 



3 f^ 

- / dcos9sin^ 9F{cos( 



3 f' 



d cos 9 cos^ 9F{cos t 



(101) 
with the function F{cos6) defined as 

/•oo 4 

F(cos 9) =. j^ dp£^ [f'{E^) + f'{E^)] , (102) 

while the transverse and longitudinal Meissner mass 
squared read 

m^ ^ e2 f^ + ^V— /' dcos9sin'9G{cos9), 
\ma mbj 4 J_i 



n 



L^ = e" I 1 I H / dcos9cos^ 9G{cosl 

\ma mbJ 2 J_i 



a_ nb\ 3e 



(103) 



(97) with the fmiction G(cos 9) defined as 



Gicose)^ ~^^y4„X/(Ep-) + /«)-! 







47r2 



CA 



47r^ \ma mb 



dp4,(^ + ^] /'«)■ (104) 



p_(3l + ^i f'fj^B 

Q Att^ \ ma mb 



In equal mass systems, we have ml^ oc px and can prove 
tya5 they are both zero 12611, which means that there are 
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no transverse Meissner effect and superfluid density. The 
reason is that the formula of the transverse Meissner mass 
squared is just the the gap equation for the LOFF mo- 
mentum, sec 4, 26]. However, for unequal mass systems, 
the gap equation seems the same as equal mass systems, 
but the formula of Meissner mass squared becomes quite 
different, and there are both transverse and longitudinal 
Meissner effects. 



VII. EXTENSION TO FINITE RANGE 
INTERACTION 



The formulae for the superfluid density and Meissner 
mass squared we derived are based on the point interac- 
tion model (0 . In this section, we show that the formula 
can be directly applied to finite range interaction sys- 
tems, if we replace the constant gap A by a momentum- 
dependent gap function A(p). 

With a finite range interaction, the Lagrangian can be 
written as 

i—a,b 

+ yd3xd3yV':(x)Vb(y)^(x,y)Vb(y)V'a(x), (i05) 

where we have assumed that the interaction is static. For 
convenience, we define the condensates 



$(x,y) = (V'h(y)V'a(x)), 
$*(x,y) = (^:(x)^,*(y)). 



(106) 



and the gap functions 



A(x,y) = V{x,y){My)M^)), 
A*(x,y) = ^(x,y)(V':(x)^,*(y)). (107) 

If the system is translational invariant with ^(x, y) = 
y(x— y), $, A and their complex conjugates depend only 
on x — y. In mean field approximation, the thermody- 
namic potential can be evaluated as 



^ = -^E/(^Trlng-i(z.„,p) 

^p)^*iq)V{p-q) (108) 



(27r)3 



(27r)3 J (27r)3 
in terms of the inverse fermion propagator 

»-'(-. p)-('lVp?~t+4)' ""' 

where V{p),^{p) and A(p) are Fourier transformation 
of y(x — y), $(x— y) and A(x — y). Since the derivation 
of the superfiuid density and Meissner mass squared de- 
pend only on the fermion propagator 5, we conclude that 
the formulae for superfluid density and Meissner mass 



squared derived in Sections and IVII are still valid in 
the finite range interaction model, if we replace the con- 
stant gap A by the momentum-dependent gap function 
A(P). 

The BP state with zero range interaction suffers neg- 
ative superfluid density, and is hence ruled out. It was 
proposed that the BP state may be stable in a finite 
range interaction model with large mass ratio 12], since 
it is the global minimum of the thermodynamic potential 
with fixed chemical potentials. For a complete study, we 
check now the superfluid density. 

For a spherically symmetric potential V{r), the gap 
function depends only on |p| and satisfies the integral 
equation 



A(g) = 



'^'Pvn . e(6g)-9(-6g) ^^^ 

V (|q - p|) — , ^^A(p). 



(27r) 



2^{^p-pif + A^{p) 



(110) 

For a given potential V and Fermi surface mismatch, we 
can solve the equation and then determine the ground 
state with the lowest thermodynamic potential. Once the 
BP solution is obtained, we can calculate the superfluid 
density and Meissner mass squared. 

For simplicity, let us concentrate on the superfiuid den- 
sity in the case with ph > Pa- The superfiuid density is 
still in the form (|70|l but the coefficients become 



as 



1 



X 



Ps=X- 



\gip^)-{X + l)vf\ 
A 

\giP2)-ix + iy,\ 



3(A- 
3(A 



1) 



dp—u 
Pi ^ 

J P 2 

dp-^Vp, 

P2 



where the function g{p) is defined as 

A(p)A'(p) 



gip) 



rrib 



P ^{^,-^if + A^{pi 



(111) 



(112) 



with A'(p) = dA{p)/dp. For the point interaction model, 
we have A'(p) — and hence g{p) = 1 which leads to 
negative superfiuid density, as we discussed above. 

Generally, the gap function peaks at p = po and drops 
down fast for p > pq. In this case, we have g{p) = 1 
and Vp = for p > P2 [13 j and the sign of the superfluid 
density depends only on a^ ■ The condition to produce a 
BP state with positive superfluid density is then 



\gipi)~iX+l)vl\> 



X 



l-3{X + l)p^'J^Upp 



'^'^u^ 



(113) 



If the slope of the gap function at p — pi is very large, 
the condition can be easily satisfied. 

For the interaction with a momentum cutoff pA^l^, the 
momentum structure of the gap function is A{p) — A for 
p < PA and A(p) = for p > pA, and we have gip) = 1 
for all p except at p — p\. Since in general case the 
positions of the zero nodes are not exactly located at the 
cutoff, pi.2 y^ PA, the situation in this model is just the 
same as in the point interaction model. 
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VIII. CONCLUDING REMARKS 

We have derived the superfluid density and the Meiss- 
ner mass squared for the fermion cooper pairing with un- 
equal masses via a standard field theory approach. For 
equal mass systems, the two variables are indeed pro- 
portional to each other, while for unequal mass systems, 
this relation breaks down. In the breached pairing states 
with zero range interaction, the superfluid density is al- 
ways negative, but the Meissner mass squared is posi- 
tive in a wide region. As a consequence, the momentum 
configuration of the LOFF pairing should be correctly es- 
tablished. We propose a proper momentum configuration 
for LOFF pairing with unequal masses and show that the 
single plane wave LOFF configuration in unequal mass 
system is physically equivalent to an anisotropic state 
with spontaneously generated superflow. These conclu- 
sions are valid only in weak coupling. Whether they are 
valid for stronger coupling, especially in the BCS-BEC 
crossover, should be examined. 

There are some problems related to the arbitrariness in 
the phase transformation induced by the U{l)^^i^U{l)^^ 
symmetry. To investigate the Goldstone mode or the 
phase fluctuation in the superfluid state, we often ne- 
glect the fluctuation of the amplitude of the order pa- 
rameter and write 0(x) = Ae^***^^^. Using the stan- 
dard phase transformation ipii^) = ipii^)^ "^6 can 
obtain the effective action for the Goldstone boson. How- 
ever, generally we can transform the fermion fields as 
ipi{x) = ■!/Ji(a:)e*'^'^^^-* with Va and Vh arbitrary constants 
satisfying the constraint i^a + f fc = 2. The low energy 
effective action for the phase field 9{x) generally reads 



'5e//[0] = -iEW-7'q^)l%o,q)|^ 



(114) 



Only when all fermionic excitations are gapped and at 
weak coupling limit, we find V and V are independent 
of Va and Vh. Hence the result of Goldstone boson ve- 
locity in 41] is safe. This problem may imply that we 
can not safely neglect the fluctuation of the amplitude of 
the order parameter in strong coupling or in the gapless 
phases. 



Another problem is the stability condition related to 
the phase fluctuation. The superfluid density ps is of- 
ten regarded as a quantity to judge the stability of BP 
state 13]. For simplicity we focus on the equal mass case. 
When Ps is negative, it directly means that the LOFF 
state has lower energy than the BP state. However, this 
is true only for the standard LOFF state with qa = qh- 
For general case with q^ = z/^q and qf, — z^^q, we should 
check the sign of Kg = d'^Vt/dq^\q=o for all possible Va 
and Vb. While ps is positive in strong coupling BEG 
region which means BP state is stable against the stan- 
dard LOFF state with qa = qb0|, there is no direct 
observation that Kq is positive for any Va and vj,, such 
as I'a = 2, j/ft = 0. If Kg becomes negative for t^Q 7^ i/(,, a 
non-standard LOFF state with q^ ^ q^, is energetically 
favored in strong coupling. 
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